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We apply a recently developed semiclassical theory of short 
peridic orbits to the stadium billiard. We give explicit expre- 
sions for the resonances of periodic orbits and for the ap- 
plication of the semiclassical Hamiltonian operator to them. 
Then, by using the 3 shortest periodic orbits and 2 more liv- 
ing in the bouncing ball region, we obtain the first 25 odd-odd 
eigenfunctions with surprising accuracy. 

PACS numbers: 05.45.Mt, 03.65. Sq, 45.05. +x 

The study of semiclassical tecniques in order to ob- 
tain quantum information of classically chaotic Hamil- 
tonian systems, has received much attention in the last 
30 years most of the developed methods are re- 

lated to the Gutzwiller's trace formula Q. This formula 
is very attractive because gives the energy spectrum of a 
bounded system in terms of periodic orbits (POs). How- 
ever, the number of POs required in the calculuation is 
enormous and increases exponentially with the Heisen- 
berg time Th = 2T:hpE {pe is the mean energy density). 

Recently, a new approach has been developed which 
uses a very little number of short POs in the chaotic re- 
gion. In order to verify the power of this new formalism, 
we applied it to the Bunimovich stadium billiard with 
radius R = I and straigth Hue 2, an ergodic system [||. 
The starting point is the construction of resonances of 
a given unstable PO. They are functions (highly local- 
ized in energy) living in a neighborhood of the PO. We 
will give the classical elements in order to obtain explicit 
expresions for the resonances and for the application of 
the semiclassical Hamiltonian operator to them. Then, 
we select a set of resonances such that its mean density 
agrees with semiclassical prescriptions. Finally, we will 
evaluate eigenfunctions and eigenvalues of the billiard by 
solving a generalized eigenvalue problem. 

Let 7 be a PO of the desymmetrized stadium billiard 
with turning points (a libration pO[|); see Fig. 0. Let x 
be the coordinate along 7 with the origin a; = at one 
of the turning points; the other being at a; = L/2, with 
L the length of 7. The transversal coordinate is y, with 
y = on 7. 

Let M{x) be a symplectic matrix describing the lin- 
earized transversal motion along 7; that is, a point with 
transversal coordinates {y,Py) at x = evolves accord- 
ing to the following rule: {y{x),Py{x)) — {y,py)M{xy . 
Then, M{x) = (-l)^(^)M(a;), with N{x) the number of 
bounces with the desymmetrized boundary while evolv- 
ing from to x, is obtained with two types of matrices: 



Mi{l) describes the evolution for a path of length I with- 
out bounces with the circle (the transversal momentum 
is measured in units of the momentum along the tra- 
jectory), and M2{9) has into account a bounce with the 
circle (9 defines the angle between the incoming trajec- 
tory and the radial direction) . ^/M^Jff) is obtained from 
M2{6) by replacing the 2 by 1. 

If the point x = (or x = L/2) is over the circle, 
we divide the contribution M2{0) given by the bounce, 
between the incoming and outgoing path. For example 
(see Fig. 0), M{L/2) = Mi{V5)^/tyWj for orbit (a), 
Mi(l + V2),/M t(0) for (b), Mi(V3)M2(V6)Mi(\/3/2) 
for (c), yM2T0)Mi(l + V5) for (d), ^AI^Mi{1 + VT0) 
for (e), and Mi(l/v^)M2(7r/4)Mi(l + 1/^2) for (f). 

By using time reversal, it is easy to see that 



M{L) 



d b 
c a 



a b 
c d 



where the explicit matrix on the right is M{L/2) and the 
other completes the orbit (from x = L /2 to x = L = 0) . 
Moreover, as the diagonal elements of M{L) are equal, 
the matrix can be written as follows H 



M{L) = i-ir 



cosh(AL) 
sinh(AL) tan((/3) 



sinh(AL) / tan((p) 
cosh(Ai) 



A = (1/L) ln{\A\ + y/A^ - 1) {A = ad + be) is the Lya- 
punov exponent in units of [length~^], tan(<y9) (7^ 0) in 
units of [length^^] defines the slope of the unstable man- 
ifold in the plane y — Py (the slope of the stable manifold 
being — tan(<y9)), and is the maximum number of con- 
jugated points along 7. Finally, being ^„ and S,s the un- 
stable and stable directions respectively, the symplectic 
matrix B transforming coordinates from the new direc- 
tions into the old ones (y and py) is 



1/a —s/a 
s a a 



B - (e. 6) = (1/V2 



M^{1) = 



1 / 
1 



and M2{e) = 



1 

-2/cos(6') 1 



with a = yj\ tan((y5)| ~ \ac/bd\^^'^, and s = sign(</j) = 
sign(acA). 

Now, it is possible to construct a family of resonances 
associated to 7. Resonances within a family are identi- 
fied by n = 0, 1, . . ., the number of excitations along the 
trajectory, and the wave number k used in the construc- 
tion depends on 7 and n through the Bohr-Sommerfeld 
quantization rule: 

Lk - {Ns + ShNh + SvNv)n - i^n/2 = 2mT. 
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V is equal to the number of bounces with the circle, Ng 
with the stadium boundary, and (Ny) with the hor- 
izontal (vertical) symmetry line. Sh = (1) for even 
(odd) symmetry on the horizontal axis and equivalently 
with Sy for the vertical axis, {v, Ns,Nh, Ny) = (1, 2, 1, 1) 
for (a), (1, 2, 2, 1) for (b), (2, 2, 2, 1) for (c), (1, 3, 3, 1) for 
(d), (1,4,4,1) for (e), and (2,2,1,1) for (f). 

Resonances are constructed with straight lines by as- 
sociating a semiclassical expression to each one. The first 
line is defined by the segment of 7 starting at xi = 0. Let 
X2 (> a^i) be the value of x such that the path reaches 
the stadium boundary while evolving along 7. The path 
going out of X2 defines the second line, and so on up to 
X = L/2. It is necessary 1 line for (a) and (b), and 2 lines 
for (c), (d), (e) and (f) (see Fig. |). 

Defining local coordinates {x^^\y'^^^) on each line such 
that x^^^ = X inside the desymmetrized billiard, the ex- 
presion for line j is [in the following expresions we are 
going to use {x,y) understanding {x^J\y^^^)] 

i^]{x,y) = fj{x,y) sm[ky'^gj{x) + kx ~ <^j{x)], (1) 

with g,{x) - [Q*{x)P,{x) + Q,{x)P*{x)]/{4 \Qj{x)\^) , 

f,{x,y) = 2(A;A)i/4exp[-fcj/V(2 \Q,{x)n/ ^ L\Q,{x)\, 
and ^j(x) = TTND{x+) + ipj + [(j)j{x)+aj{x)]/2. Nd{x+) 
is the number of bounces up to Xj (including the 
bounce at Xj) satisfying Dirichlet boundary conditions; 
Nd{xi) = 0. Moreover, 

with 

Xj' means (for j > 2) that M is evaluated after the bounce 

with the boundary at xj [M{x^) = M2(7r/6)Mi(y3/2) 
for orbit (c), Mi(\/5/2) for (d), Mi(2\/TU/3) fo r (e), an d 
M2(7r/4)Afi(l + 1/V2) for (f)]. M{xt) = y/th{0) if 
xi lies on the circle; otherwise M{xi) — 1. 4'j{x) = 
a.rg[Qj{x)] — arg[Qj(a;j)] (arg takes the argument of 
a complex number in the range [0, 27r)). aj{x) = 
27r sign(a; — Xj) if nj{x) ^ nj{xj) and {x — Xj)4>j{x) < 0; 
otherwise aj{x) = 0, where 

{1 if X — > max(— Cj/cj, —bj/dj) 
— 1 if X — Xj < min(— a_,/cj, —bj/dj) 
otherwise. 

If Cj =0 or dj =0, Xj is replaced by any other point on 
line J, inside the desymmetrized billiard. (j)j{x) A- aj{x) 
defines the angle swept by Qjlx'^'^^) in a continuous way. 
Finally, (pj = ipj-i + [0j_i(xj) + aj^i{xj)]/2 for j > 2. 
The value of ipi depends on the starting point and the 
symmetry, and all the possibilities are considered in Fig. 



|T|. ipi is equal to — s/i7r/2 for (a), for (b), {sh — l)7''/2 
for (c), (s/i + St, — l)7r/2 for (d), —Sv1t/2 for (e), and 
(s„ - l)7r/2 for (f). 

The transformation from local coordinates (x'--'-' , y^^^ ) 
on line j to coordinates {X, Y) (horizontal and vertical 
directions in the plane respectively) is obtained through a 
simple transformation. If {Xj,Yj) are the coordinates of 
the point Xj , and aj the angle of line j with the horizontal 
direction, (x^-'-' — Xj,y^^^) — Gj{X,Y) is given by 

G,{X, Y) = {X-X„Y- Y,) ( '"^f^l - '''l^''^^ ) . 
■'^ ' ^ \^ sm(Q!j) cos(aj) J 

Finally, the family of resonances i^^{X, Y) is constructed 
with all the lines including symmetries (see Fig. H) 

^7 = E E E ^' [{xj , 0) + G, {siX, s,Y)] . (2) 
j 1=1 1=1 

= i-iy+' and si = {-iy+\ h, = [^.,1 +<5,,2(1 - 2s,,)] 
and vi — [Si^i + (5;, 2(1 ~ 2s^)]. rrih and rriy depend on 
j and are specified as follows: m/i = 1 (2) if the line is 
(is not) symmetric with respect to the horizontal axis, 
and equivalently with my for the vertical axis; however, 
rufi — 2 and my = 1 if the line goes through the origin. 

We define H = — and E = (remember that k 
depends on 7 and n). Then, the semiclassical approx- 
imation for {H — E)^p^{X,Y) is obtained directly from 
Eq. ^ having into account the following semiclassical 
prescription [|| 

iH-E)^,{x,y) f,{x,y) sm[ky^g,ix) + A{x)l (3) 

with fjix,y) = Xk (2fcyVlQiWP - ^)fj{x,y), and 
A(x) = kx - $j(x) + tt/2 ~ 2aig[Qj{x)]. That is, the 
action H — E on tpj excites the transversal direction 
with two excitations. Using these expresions it is pos- 
sible to obtain matrix elements by direct integration on 
the domain (the quarter of billiard in this case) |l^ . 

In the semiclassical limit (n 00), the following di- 
agonal matrix elements are obtained explicitly | p[| (us- 
ing Dirac's notation): i) (7I7) 1 , ii) E = 

(7|i7|7)/(7|7) ^ , and m) a' = {j\H^\j) / - 
2 

E 2A^fc^, with a the dispersion of ip-y. On the other 
hand, as the operator H is not exactly Hermitian, it is de- 
fined a symmetrized interaction between POs as follows: 
{6\H\j)^i{S\Hj) + {j\HSy)/2. 

Now we will select the set of resonances defining an 
appropiate basis. From now on, we will restrict to odd 
symmetry on the horizontal and vertical directions. For 
the selection of bouncing ball functions we use Tan- 
ner's prescription p^ . That is, for a given number 
M (M = 1,2,...) of vertical excitations, there are 
N — [\/2M + 1 + 5/4] horizontal excitations in the range 
M < k/ir < AI + 1. This translates into a number n of 
longitudinal excitations (along the selected orbit) that is 
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specified by n = {M—1)K+N—1, where K is the number 
of segments of the trajectory inside the desymmetrized 
billiard [for instance, K = 2 ioi (b)]. Then, for the con- 
struction of the function (M, TV) we select the trajectory 
(among (b), (d) and (e)) which gives the resonance with 
smallest dispersion a. 

The number of bouncing ball wavefunctions is only a 
fraction of the mean number of states specified by the 
Weyl's law. Then, using the shortest periodic orbits (a), 
(b) and (c), we put as many resonances as necessary to 
obtain the required number. As the period of the three 
POs is comparable, we select first orbit (c) because the 
associated resonances have smallest dispersion. Figure ^ 
clarifies the situation. The first column shows the spec- 
trum of bouncing ball resonances. Over each line appears 
the label corresponding to Fig. §(a), M, N, the orbit 
used, and n. Second column shows the spectrum of res- 
onances constructed with orbit (c), and over each line 
appears the label and n. The same is applicable to the 
third and fourth columns with orbits (a) and (b) respec- 
tively. In this way, the density of resonances agrees with 
the semiclassical mean density; however for k > 15 more 
orbits are required. Note that orbit (b) is used for the 
construction of resonances living in the bouncing ball and 
chaotic regions. This is because the bouncing ball region 
decreases as l/Vk. Orbit (f) is not considered because 
the associated resonances do not satisfy boundary condi- 
tions with sufficient accuracy for low energies. 

Semiclassical eigenfunctions are constructed with this 
set of wavefunctions. Of course, only a limited number 
of them are required for a particular eigenfunction. State 
21 (see Fig. ^(b)) needs resonances from 18 to 23 (Fig. 
^(a)). The other states use less resonances. Suppose that 
at fco there is an eigenstate of the system; it is constructed 
with all resonances satisfying \k — ko\ < 0.8 Despite 
this is a semiclassical criterium, it works in general at low 
energies too. It says that the number Nr of resonances 
contributing to each eigenstate increases as follows: Nr ~ 
0.5 fco- 

Figure ^(a) shows linear density plots of the 27 se- 
lected resonances arranged by energy. Numbers below 
each plot are the label (left), the squared root of the mean 
energy, and the dispersion a in units of the mean energy 
spacing. Using this basis of functions, we evaluate the 
overlaps and the Hamiltonian matrix elements using the 
semiclassical prescription given in (^. Then, by solving 
a generalized eigenvalue problem, the semiclassical set of 
eigenfunctions shown in Fig. ^(b) is obtained. Num- 
bers below are the label, the semiclassical wave number 
and the overlap with the exact solutions, which arc dis- 
played in Fig. ^(c). The mean standard deviation of 
the semiclassical eigenvalues is a fraction (0.06) of the 
mean level spacing in accordance with the theory. Re- 
cent results obtained for the hyperbola billiard ||l^ show 
a standard deviation of 0.047 (in units of the mean level 
spacing) for the first 24 even eigenvalues. Those results 



were calculated using trace-formula-type techniques, in- 
volving 38 131 periodic orbits (193 695 pseudo-orbits), in 
contrast with the only 5 used in the present work. On the 
other hand, overlaps of the semiclassical eigenfunctions 
with the exact ones are surprisingly good (see Fig ^(b)). 
Moreover, to our knowledge, this is the first semiclassical 
evaluation of a set of eigenfunctions in a chaotic system 
using periodic orbits. 

In conclusion we have applied successfully the theory 
of short POs to the stadium billiard. This shows that the 
classical information contained in short POs is sufficient 
for obtaining the stationary states of a bounded chaotic 
Hamiltonian system. 
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03-00050-01015, SECyT-ECOS and CONICET PIP 
0420/98. E.V. acknowledges the hospitality of the 
LPTMS, Orsay, where part of this work was done. 

FIG. 1. Set of periodic orbits of the desymmetrized sta- 
dium billiard used (with the exception of orbit (f)) for the 
construction of resonances. (1) and (2) labels the correspond- 
ing straight line. 

FIG. 2. Set of lines including symmetries used for the con- 
struction of resonances associated to orbit (d). The different 
sets of coordinates used are indicated. 

FIG. 3. Spectrum of resonances used for the evaluation 
of semiclassical eigenfunctions. The meaning of each column 
is as follows: the first number over the lines is the label. In 
the first column, showing the spectrum of bouncing ball reso- 
nances, M, N and n are indicated. Letters correspond to the 
orbit from which they were constructed. In the other three 
columns only the label and the number n is over the lines; the 
corresponding orbit, the same for the entire column, is at the 
bottom of the figure. 

FIG. 4. Linear density plots of: (a) the selected reso- 
nances. The numbers below each plot are (from left to right) 
the label, the squared root of the mean energy, and the dis- 
persion a in units of the mean energy spacing, (b) the semi- 
classical eigenfunctions. Here, after the label, the semiclas- 
sical wave number and the overlap with the corresponding 
exact eigenfunction are displayed, (c) exact eigenfunctions 
with their wavenumbers. 
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(27) 4,4 (d) 12 



(23) 4,3 (d) 11 


(21) 


4,2 (d) 10 


(18) 4,1 (e) 12 


(14) 


3,3 (d) 8 


(12) 


3,2 (d) 7 


(10) 


3,1 (e) 8 


(8) 


2,3 (b) 4 


(6) 


2,2 (b) 3 


(4) 


2,1 (d) 3 


(2) 


1,2 (b) 1 


(1) 


1,1 (b)0 



M,N n 



(25)9 
(20) 8 
(16)7 
(13)6 
(9) 5 
(7) 4 
(5) 3 
(3) 2 



(24) 8 
(19)7 
(15)6 
(11)5 

n 



(26)9 
(22) 8 
(17)7 

n 



n 

(c) (a) (b) 



